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■s^j- ' Let g be the Lie algebra of a compact Lie group and let 9 be any automorphism of q. Let t 

£<\ , denote the fixed point subalgebra g e . In this paper we present LiE programs that, for any finite 

dimensional complex representation n of q, give the explicit branching 7r|{ of 7r on t. Cases of 
special interest include the cases where 9 has order 2 (corresponding to compact riemannian 
symmetric spaces G/K), where 9 has order 3 (corresponding to compact nearly-kaehler homo- 
geneous spaces G/K), where 9 has order 5 (which include the fascinating 5-symmetric space 
■ E^/A^Ai), and the cases where 6 is the centralizer of a toral subalgebra of g. 
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1 Introduction 

There are many situations where one wants to see the explicit branching of a particular represen- 
tation from the Lie algebra q of a compact Lie group to a Lie subalgebra 6. In many cases the 
situation corresponds to a compact homogeneous space G/K of some geometric interest, such as 
the cases where G/K is a riemannian symmetric space, a nearly-kaehler manifold, or the compact 
group realization of a complex flag manifold. Most cases of geometric interest have the interesting 
property that 6 is the fixed point set of an automorphism 6 of g. In essentially all cases one can 
compute the branching by hand, but the time and effort involved may be extreme. This situation 
is greatly ameliorated by use of the public domain computer program LiE [9]. In this paper we 
produce the LiE routines that carry out the branching of representations from q to t explicitly 
' when 6 is the fixed point set of an automorphism 6 of q. 

a 

One might expect the built-in branch routine of LiE to do the job for us without any additional 
programming. The problem is that LiE mixes up the order of simple roots, making iteration of 
branching very difficult and causing serious problems for identifying the restriction in cases where 
there is a symmetry of the Dynkin diagram of t. Worse, on each summand of the restricted 
representation it renormalizes the restriction to the center of t in a complicated manner, and that 
causes even more serious problems in geometric and analytic applications where negativity is needed 
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and is controlled by restriction to the center of t. Our LiE routines specifically address and solve 
those problems. 

We developed many of these LiE routines for use in our work [S] on the range of the double 
fibration transform [6l Chapter 14], where we need explicit information on branching from the Levi 
component of a parabolic subgroup to its intersection with a maximal compact subgroup. These 
LiE programs are based on structural information on Lie algebras and automorphisms to be found 
in [2], [IDl, [IT], m and 0. 

We necessarily start out by describing use of the LiE program and how its use varies with the 
properties of (q,0). Thus in Section [2] we indicate root orderings and their role in computing LiE's 
"restriction matrix" . Then in Section I2AI we reduce questions of branching to the cases where g 
is simple and t is a maximal ^-stable subalgebra of g, where there are three essentially different 
situations. The case where g is simple and rankt < rankg is described in Section I2BI It relies 
on information from [10], [11] and [Tj- The case where g is simple, rankt = rankg and t is not 
semisimple, is the subject of Section I2CI It relies on information from [2], [7], and the standard 
structure theory of parabolic subgroups. Then the case where g is simple, rank 6 = rankg and 6 is 
semisimple, is indicated in Section I2D1 This is the most delicate case, and it depends on methods 
from [2], [7] and [8]. 

In Section [3] we list all cases where g is simple, 6 is ^-maximal and rank A: < rankg. For each of 
them we describe how to find the restriction matrix and we give the listing of a LiE program that 
computes branching from g to t. The programs are (|3.ip . (|3.2p . (|3.3p . (|3.4p . (|3,5p . (|3.6p . (|3.7[) and 
()3.8p . In all but two of these, 6 2 = 1 so G/K is a riemannian symmetric space, and in those two 
we have 6 3 = 1. 

In Section H] we discuss the LiE programs for the cases where g is simple, rank 6 = rankg and 
t is not semisimple. Those essentially are the cases where g is simple and t is the centralizer of a 
toral subalgebra, where the LiE programs are described in Section I2CI 

Section [5] gives the LiE branching programs for the cases where g is simple, 9 2 = 1 and rank t = 
rankg. The programs (|5.ip . (|5.2p and (|5.3p apply when g is classical. There one has no surprises 
on the root orders, but when g is exceptional the LiE program scrambles the root order going from 
to t In ([S3D , (pT5j) . (joTJj) . (pU2|) . pTT3|) and (|5T4"j) this is fairly straightforward, as there is not 
much flexibility for the location of t inside 0. However, in applications [5] we need to keep track 
of the various simple roots, and we must deal with the fact that there are three combinatorially 
distinct ^1^5 's in Eq and two essentially distinct A\D§s in Ey. This results in more programs 
than one might expect, specifically in ([521), d5HD , <^Mi dSHOD and flSTTT]) . 

Section [6] completes the results of Section [5j providing the LiE routines for the seven remaining 
cases, those where is simple, 6 3 = 1 or 9 5 = 1, and rankfi = rank0. These routines are (16. ip . 
(|6.2p . (|6.3p . (|6.4p (|6.5p . (|6.6p and (|6.7|) . There, as in the exceptional group cases of Section we 
label the simple roots of t to minimize any departure from the root ordering of g. 

As indicated in Section I2A[ this completes the analysis of branching of finite dimensional irre- 
ducible representations from the Lie algebra g of a compact Lie group to the fixed point set t of 
any automorphism 8 of g. 
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2 Restriction Matrices and Branching in LiE 



All our LiE routines are given by files with names of the form branch_X_Y. lie where X is the 
LiE designation of the type of q, e.g. E6, and Y is the LiE designation of the type of t, e.g. F4. 
They are called within the LiE program by first reading in the file, (> read branch_X_Y. lie) and 
then giving the command (> branch_X_Y(v) ) where v = [v\, . . . ,v n ] represents the highest weight 
^2 Vi^i of an irreducible representation ir of q to be branched on I. Here the £j are the fundamental 
simple highest weights. Note that this depends on the ordering of the simple roots tp^. LiE uses 
(and therefore we use) Bourbaki order [3j, given as follows on the Dynkin diagrams. 



02 



tpi 



(type A e , £ ^ 1) 



o o- 



01 02 



o O- 



01 02 



o o- 



01 02 



- o ) o 



ipt-\ tpt 




(type B e , £^2) 
(type C t ,lZ3) 

(type B t , £ ^ 4) 



q \ 3 3 



o o ) o o 

Ipl 1p2 1p3 1p4 



(type G 3 ) 
(type F 4 ) 



Ipl 1p3 4>4 i>5 4>6 

o o o o o 



6 ^2 

Ipl 1p3 "04 i>5 -06 07 

o o o o o o 



6 02 

01 03 04 05 06 07 

o o o o o o — 



6 02 



08 
^3 



(type E 6 ) 



(type E 7 ) 



(type E 8 ) 



where, if there are two root lengths, the arrow points from the long roots to the short roots. 

If 6 is a subalgebra of q then the LiE program computes branching of representations by use of a 
"restriction matrix" . This is the matrix whose rows are the restrictions, from a Cartan subalgebra t 
of g to a Cartan subalgebra s C t of t , of the fundamental simple weights of q as linear combinations 
of the fundamental simple weights of t. Obviously this depends on the relation between our choices 
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of simple root systems for q and 6. 



2A Reduction to the cases where g is simple and £ is maximal. 

We start with the Lie algebra of a compact connected Lie group G and an automorphism 9 of 0. 
The fixed point algebra is t = q 9 , and K is the corresponding analytic subgroup of G. We start 
also with an irreducible finite dimensional representation ir of g. We want to describe tt\(_ explicitly. 

We indicate how to reduce our branching questions to the case where q is simple and t = q 9 
is maximal among the ^-stable subalgebras of q. That done, we have three essentially different 
possibilities. The methods appropriate to those three situations are addressed in Sections I2B} I2CI 
and I2DI below, and carried out completely in the remainder of this paper. 

Our branching procedures all use the LiE program. We give listings of the relevant LiE routines, 
and when the programming aspects are not so obvious we give an exposition of the mathematics 
behind our branching routines. 

Write q = q' © 3 where q' is semisimple and 3 is the center of 0. Each summand is 0— stable, 
so 6 = (t n q') ©(603). Also 7r = 7r' M x, exterior tensor product, where tt' represents t n g' and 
X is a 1 -dimensional representation of 3. Now 7r|j = (7r'|tng') Kl (xknj) an d evaluation of the latter 
factor is just restriction of a linear functional to a linear subspace. Thus we need only worry about 
computing vr'|f ng /. That is the first reduction: it suffices to consider the case where is semisimple. 

Decompose g as a direct sum of simple ideals. Then 9 gives a permutation on that set of 
ideals, and as such it is a product of disjoint cycles. In other words, we have a decomposition 
g = f)i © • • • © f)r where 9 preserves each f)j and induces a cyclic permutation on its simple direct 
summands. Now 6 = q 9 = t)f © • • • © t)f. That is the second reduction: it suffices to consider the 
case where 9 induces a cyclic permutation on the simple ideals of g. 

Now we have reduced to the case g = gi © • • • © g m where the g^ are simple, #(gi-i) = 0i for 
1 < i ^ m, and 9(g m ) = gi. We interpret the 9 : gj_i = g^ as identifications. That done, 

(2.1) g = 01 © • • • 01 (m summands) where 0(£i, ...,£ m ) = (7(£m),£i, • • • ,£m-l) for & G gi. 
Here 7 is an automorphism on gi. Now we have 



Thus 6 = q 9 = (q1 © • • • © qI) 9 where there are m summands gj. Denote 61 = 0i- From (12. ip and 
(j2T2|) we have 



Now it suffices to consider the case where is simple. 

We address the programming aspects. Suppose that we are given an irreducible representation tt 
ofg = gi©---©gi (m summands) . Then tt is the exterior tensor product 7Ti M ■ ■ ■ M ir m of irreducible 
representation 7Tj of 01. In view of (j2.3fl . is the interior tensor product of the restrictions of the 
7Tj to the 61 = gj. We can do this in two stages. First we compute the restrictions vrj^, which 
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only involves cases where we branch from a simple Lie algebra, and then we decompose the tensor 
product. In the latter setting we have reduced to the case where 7 = 1 but gj may no longer be 
simple. Still, g d decomposes under the action of 9 in the setting of a cycle of simple ideals. This is 
the third reduction: the branching problem is reduced to the case where = 0i © ■ ■ • © 0i, sum of 
m simple ideals, and 9 acts by #(£1, . . . , £ m ) = (£ m , £i> • • • > Cm-l)- I n this case t = ® • • • <8> vr m 
and 6 is the diagonal diaggi in g. 

We have reduced the case of branching from non-simple g to two parts: branching from simple 
proper subalgebras of g and decomposing tensor products of irreducible representations of g. The 
latter is done in LiE as follows. Let v be a matrix of m rows, each row v[i] a vector of length 
equal to the rank n of 0i, where the row v [i] = [v [i , 1] , ... , v [i ,n] ] describes the highest 
weight Aj = Y2j v [hj]£j °f n i i R terms of the fundamental simple weights £j. If m = 2 and the 
default is set to the Cartan type of gi then we can use LiE's built-in function 

tensor (v [2] , v [1] ) 

for the tensor product decomposition. If m > 2 we do this recursively, but we must first convert 
the v[i] to polynomials in the LiE sense, 

w = null(m,n); for i = i to m do w[i] = tensor (v [i] , null (n)) od 

and then we can issue the LiE command 

w = w[l] ; for i = 2 to m do w = tensor (w[i] ,w) od 

Here is a general LiE routine to systematize this. It is called in LiE by branch_diag(v,g) where g 
is the Lie type of a simple Lie algebra such as A3, C7, G2, F4 or E8, and where v is a matrix of non- 
negative integers whose rows have length rank g representing highest weights of the representations 
of g to be tensored together. 

# file branch_diag . lie # 

# usage: branch_diag(v,g) where g is a simple Lie algebra type # 

# (An, . . . , E8) and v is a matrix of rank g columns, whose rows # 

# specify the highest weights of reps $\pi_i$ of g; It returns # 
, # the (interior) tensor product of the $\pi_i$. # 

branch_diag(mat v; grp g) = setdef ault (g) ; 

loc u = tensor (null (Lie_rank) , null (Lie_rank)) ; 

for r row(v) do u = tensor(u,tensor(r,null(Lie_rank))) od; 

print ("the branching from product of "+n_rows(v)+" copies of " 

+Lie_group(Lie_code [1] ,Lie_code [2] )+" to the diagonal is"); u 



2B Case g simple and rankt < rankg. 

Suppose first that is simple and rankt < rankg. Choose respective Cartan subalgebras s C t. 
Then there is a simple root system ^ = {ipi, . . . , tp n } for (g, t) such that the restrictions ipi\ s , ■ ■ ■ ipn\s 
form a simple root system $ = {(pi, . . . , (f r } for t. See |10j . In that case we have a root restriction 
matrix res_rt whose j th row is given by res_rt[j] = [mji,...mj r ] where ipj\s = '^k m j,k l Pk- 
LiE however requires the corresponding restriction matrix of fundamental simple weights, and can 
compute it from res_rt as 

(2.5) res_wt = i_Cartan(g)*res_rt*Cartan(t)/det_Cartan(g) 
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where i_Cartan(g) /det_Cartan(g) is the inverse of the Cartan matrix of g using ^ and Cartan(6) 
is the Cartan matrix of t using 

Here is an example. Let g = su(7) and 6 = 50 (7). The Dynkin diagram of t is obtained by 
folding that of g, 
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In other words, the simple root restrictions are tpi i— ► (pi, ip2 l— *• <J>i-> V>3 



and ^6 l— > Pi- Thus res_rt is 



10 
10 
1 
1 
10 
10 



. Now (|2.5p gives res_wt 



10 
10 
2 
2 
10 
Ll 0. 



> ip 3 , -04 ¥>3, 0>5 ^ ^3 

If the LiE default group 



is set to A6 for su(7) then branching of the adjoint representation of su(7) on so (7) is given by 
branch([l, 0,0, 0,0,1] ,B3,res_wt) , resulting in 1X[0,1,0] +1X[2,0,0]. 



2C Case g simple, rank t = rankg and t is not semisimple. 



Suppose that g is simple and 6 is of equal rank but is not semisimple. Recall that t is ^-maximal 
in the sense that it is maximal among the ^-stable proper subalgebras of g. It follows that t is 
the centralizer of its center, so it is a compact real form of the reductive (Levi) component of a 
parabolic subalgebra of gc. We remark that the centralizer of a toral subalgebra e C t of g is always 
the fixed point of an automorphism 9 £ Ad(exp(e)), for example 9 = Ad(t) where the powers of 
t form a dense subgroup of the torus exp(e). Now g has a simple root system ^ = {ipi, . . . ,ip n } 
such that some subset C \£ is a simple root system for t. We use the notation of Baston & 
Eastwood [llj to indicate these Levi components, i.e. to indicate these centralizers in g of subtori 
of its Cartan subalgebra. Thus if ip S f \ $ we replace the circle o for if) by a cross x. We 
refer to this as the diagram of the corresponding parabolic subalgebra q$ of gc, the corresponding 
parabolic subgroup Q$ of Gc, and our algebra t = q$ n g which is a compact real form of the 
Levi component of q$. For example, the parabolic subalgebra q<j> of sl(n + 1;C) that corresponds 
to the complex projective space P n (C) = SL(n + 1;C)/Q$ is given by ^ = {ipi, . . . ,ip n } and 

= {"01, . . . , -0n-i}) so it has diagram ° ° — ■ •• — x , and the parabolic subalgebra 

for the Grassmannian of lines in hyperplanes in C n+1 is given by $ = {ip2, ■ ■ ■ > VVi-l} an d has 

diagram x — ■ ■• — ° x . These correspond to the cases t = u(n) C su(n + 1) = g and 

t = {i6 u(l) © u(n - 1) © u(l) | trace x = 0} C su(n + 1) = g. 

Suppose that consists of all but one element 7 = tpi of in other words that q<j> is a maximal 
parabolic subalgebra of gc- That is the case where there is only one x on the diagram of q$. Then 
there is a simple LiE routine (from the LiE manual |9j) that describes branching of representations 
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from g to t = q$ n g: 

# file Levi_branch. lie # 

Levi_mat (int i) = f undam(id(Lie_rank) - i) 
Levi_type (int i) = Cartan_type (Levi_mat (i) ) 
, . Levi_diagram(int i) = diagram(Levi_type (i) ) 

Levi_res_mat (int i) = res_mat (Levi_mat (i) ) 
Levi_branch(vec v; int i) = loc m = Levi_mat (i) ; 
r = res_mat (m) ; 
branch(v, Cartan_type (m) , r) 

We use it, say with q = E7 and 6 = EqTi, as follows. Do read (Levi_branch. lie) , then 
setdef ault (E7) , then diagram in LiE to see that 7 = ^7) do v = [v_l, v_2, v_3, v_4, v_5, 
v_6, v_7] for the highest weight Yl v i£i °f ^1 anc ^ compute the restriction by Levi_branch(v,7). 
The result is a sum of vectors with multiplicities, e.g. IX [0 , , , , , 2 , -6] +2X [0,0,0,0,0,2,-4] 
+4X [0,0, 0,0, 0,2, -2] + ... where the last entries (-6, -4, -2) refer to the central torus. For 
the meaning of the others do Levi_diagram(7) in order to compare the root orderings (in the LiE 
program) for q and 6. 

Suppose next that $ consists of all but two elements ipi and ipj of VP, in other words that there 
are two x's on the diagram of q$. We modify the routine (|2.6p to accommodate this. Here it is 
important that i > j so that we remove rows i and j from a matrix by removing the i th and then 
the j th of that. 

# file Levi_branch2 . lie # 

Levi_mat (int i, j) = f undam( (id(Lie_rank) - i) - j) 
Levi_type (int i, j) = Cartan_type (Levi_mat (i , j ) ) 
,^ _n Levi_diagram(int i, j) = diagram(Levi_type (i , j ) ) 

Levi_res_mat (int i, j) = res_mat (Levi_mat (i , j)) 
Levi_branch2(vec v; int i, j) = loc m = Levi_mat(i, j) ; 
r = res_mat (m) ; 
branch(v, Cartan_type (m) , r) 

Similarly if $ consists of all but three elements ipi, t/jj and ipk of \E', i > j > k, 

# file Levi_branch3 . lie # 

Levi_mat (int i, j, k) = f undam( ( (id(Lie_rank) - i) - j) - k) 
Levi_type (int i, j, k) = Cartan_type (Levi_mat (i , j ,k)) 
, . Levi_diagram(int i, j, k) = diagram(Levi_type(i, j ,k)) 

Levi_res_mat (int i, j, k) = res_mat (Levi_mat (i , j ,k)) 
Levi_branch3(vec v; int i, j, k) = loc m = Levi_mat(i, j ,k) ; 
r = res_mat (m) ; 
branch(v, Cartan_type (m) , r) 

At this point the pattern is clear. For example, try 

read Levi_branch3 . lie 
setdef ault (E8) 

Levi_branch3 ([1,0,0,0,0,0,0,11,8,6,4) 

The first five entries in each of the resulting 8-tuples gives the branching on the semisimple part 
A2A1A1A1 of 6 = q^\{,/, gj ^ 6j ^ 4 }, but with some roots permuted. To see the permutation look at the 
restriction matrix Levi_res_mat (int 8, 6, 4) and remove rows 8, 6 and 4, and remove the last 
three columns. In LiE this can be implemented as 
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*(((*(((Levi_res_mat(8, 6, 4) - 8) - 6) - 4) -8) -7) -6) 

In this way the Levi_branch LiE routines give the restriction to the semisimple part of 6. 

Of course, these routines also give the action of the center of 6 on each irreducible summand 
but, unfortunately, this is implemented in a rather ad hoc fashion. We now explain how to specify 
the central action in a more systematic and useful manner. At the same time, we avoid having 
to deal with the permutations introduced by the programs Levi_branch, as above. The problems 
with these programs can be illustrated with the following simple examples. With setdef ault (F4) 
in place we have LiE calculate the following matrices. 



(2.9) 



i_Cartan 


Levi_res_mat (3) 


Levi_res_mat (4) 


Levi_res_mat (4 , 3) 


[[2,3,4,2] 


[[0,1,0,4] 


[[1,0,0,2] 


[[0,1,2,0] 


, [3,6,8,4] 


, [1,0,0,8] 


, [0,1,0,4] 


,[1,0,4,0] 


, [2,4,6,3] 


, [0,0,0,6] 


, [0,0,1,3] 


, [0,0,3,0] 


, [1,2,3,2] 


, [0,0,1,3] 


, [0,0,0,2] 


, [0,0,0,1] 


] 


] 


] 


] 



In this particular case, the matrices Levi_res_mat (i) are easy to understand. The first three 
columns specify a permutation of the uncrossed nodes and the last column is the i th column of the 
inverse Cartan matrix i_Cartan. It is easy to check that the element of the Cartan subalgebra tC j 
defined by the i th column of the inverse Cartan matrix with respect to the basis of fundamental 
weights is in the center of the corresponding Levi subalgebra t. (Indeed, this is minus the so-called 
'grading element' of the corresponding maximal parabolic subalgebra [4].) Thus, the restriction 
matrix specifies a permutation of the uncrossed nodes and a particular element of the center. Here 
is the branching of the adjoint representation given by Levi_branch( [1,0,0,0] ,3). 

[0, 1,0, 4] © [0, 0, 1, 3] © [1, 0, 2, 2] © [0, 1, 1, 1] © ([0, 0, 2, 0] © [1, 1, 0, 0] © [0, 0, 0, 0]) 

© [1, 0, 1, -1] © [0, 1, 2, -2] © [0, 0, 1, -3] © [1, 0, 0, -4], 

(where the ordering is given by the value of the grading element from —4 to 4). In other words, 
the Lie algebra q = F4 decomposes as 

= 0-4© 0-3 © 0-2 © 0-1 © 00 © 01 © 02 © 03 © 04 = [0, 1, 0, 4] © [0, 0, 1, 3] © • • • © [1, 0, 0, -4] 

(and this is exactly the realization of as the |4|-graded Lie algebra corresponding to the parabolic 
subalgebra q$ = 0o © • • • © 04 as in [U Theorem 3.2.1]). 

The restriction matrix Levi_res_mat (4,3) in f|2.9[) is more difficult to understand. Certainly, 
we could use 

[[0,1,4,2] 
, [1,0,8,4] 
, [0,0,6,3] 
, [0,0,3,2] 
] 

as a more easily understandable restriction matrix. It is obtained by using the j th and i th columns 
of i_Cartan to replace the last two columns of r, a change that is easily implemented in LiE by 
adding 
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for k = 1 to Lie_rank do r [k,Lie_rank-l] = i_Cartan[k, j] od; 
for k = 1 to Lie_rank do r [k,Lie_rank] = i_Cartan [k, i] od; 

as the penultimate two lines of Levi_branch2 . lie. In comparison with (|2,9p . the last two columns 
of Levi_res_mat (4,3) are some linear combination of the appropriate columns of the inverse 
Cartan matrix. Moreover, the case g = F± is deceptively simple because its Cartan matrix has 
unit determinant. In general, because LiE is restricted to integer arithmetic, it is only reasonable 
to use the appropriate columns from i_Cartan, as above. In particular, the grading element will 
not be simply minus the sum of these columns but, in addition, one must divide by det_Cartan. 
Although the grading element takes on integral values on the adjoint representation (from —A; to 
k where g is |fc|-graded by the parabolic subalgebra q$), for a general irreducible representation 
its values will be rational with det_Cartan as denominator. In any case, the raw instructions 
Levi_res_mat (i , j) and Levi_res_mat (i , j ,k) produce rather bizarre changes of basis from the 
more natural normalization provided by the inverse Cartan matrix and even Levi_res_mat (i) is 
better modified by 



for j = 1 to Lie_rank do r [j ,Lie_rank] = i_Cartan [j , i] od; 



to avoid spurious factors. 

For many purposes, however, it is better to write all weights as linear combinations of the 
fundamental weights of (g, t) and, following [lj, attach the resulting coefficients to the corresponding 
nodes of the Dynkin diagram. In our example, the adjoint representation LJM— 2 decomposes as 



(2.10) 



0\-l 2 



1 Ov-l 1 




11-2 



o — o 



The conversion between these two conventions is the definition of the restriction matrix. Therefore, 
no matter what restriction matrix is used, to convert back to the conventions of pQ, one simply 
needs to invert the restriction matrix and apply this inverse matrix by right multiplication to each 
term obtained from branch (v , Cartan_type (m) , r) . Since LiE allows only integer multiplication, 
inverting a matrix with integer entries requires some care. In the following program, the restriction 
matrix r is inverted by first extracting from it a permutation matrix p, noting that permutation 
matrices are orthogonal, and forming (*p)*r. The result necessarily has the form 



1 

1 







* 

* 

: * 

1 * 
* 
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which is inverted by a dint of an explicit formula. 



# file Levi_branch_improved. lie # 

Levi_mat (int i) = f undam(id(Lie_rank) - i) 

Levi_type (int i) = Cartan_type(Levi_mat(i)) 

Levi_diagram(int i) = diagram (Levi_type ( i) ) 

Levi_res_mat (int i) = res_mat (Levi_mat (i) ) 

div(pol p;int k) = loc 1 = OX(null(n_vars(p))) ; 

for i=l to length(p) do 1 = 1+coef (p, i)X(expon(p, i)/k) od; 1 

Levi_branch(vec v; int i) = loc m = Levi_mat (i) ; 

loc r = res_mat (m) ; loc p = r; 

for j = 1 to Lie_rank do p[j ,Lie_rank] = od; p[i,Lie_rank] = 1; 
loc q = (*p)*r; loc det_q = q[Lie_rank,Lie_rank] ; 
loc qq = q; 

for j = 1 to Lie_rank do qq[j ,Lie_rank] = -q[j ,Lie_rank] od; 

qq[Lie_rank,Lie_rank] = 1; 

loc s = null(Lie_rank,Lie_rank) ; 

for j = 1 to Lie_rank do s[j,j] = det_q od; 

s [Lie_rank,Lie_rank] = 1; loc i_q = qq*s; 

loc b = branch(v, Cartan_type (m) , r) ; div(b*i_q*(*p) ,det_q) 



The program is used as before but the result is 
of pQ. For example, 

read Levi_branch_improved . lie 
setdef ault(F4) 
Levi_branch( [1,0,0,0] ,3) 



expressed using the diagrammatic conventions 



gives 

1X[0, 0,-1,1] +1X[0,0,-1,2] +1X[0,0, 0,0] +1X[0,0, 0,1] + 

1X[0, 1,-2,0] +1X[0, 1,-2,1] +1X[0, 1,-2,2] +1X [1 , , -2 , 2] + 

IX [1, 0,-1,1] +1X[1,0, 0,0] +1X[1, 1,-2,0] 

as in (|2.10j) (but devoid of the convenient ordering there. The ordering of (I2,10p is essential when 
the action of the full parabolic q$ is considered rather than just its Levi factor. Representations of 
q$ are generally filtered. For example, the tail of (|2.10p . 

01 -2 1 10 -2 2 00 -1 1 01 -2 
O 1 1 \ * O © O 1) \ * O © O 1 1 \ * O © O 1 1 \ * o 

is interpreted in [U p. 135] as inducing the cotangent bundle on the corresponding generalized flag 
manifold.). 
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Levi_branch2 . lie is similarly improved 

# file Levi_branch2_improved . lie # 

Levi_mat (int i, j) = f undam( (id(Lie_rank) - i) - j) 
Levi_type (int i, j) = Cartan_type (Levi_mat (i , j ) ) 
Levi_diagram(int i, j) = diagram (Levi_type (i , j ) ) 
Levi_res_mat (int i, j) = res_mat (Levi_mat (i , j)) 
div(pol p;int k) = loc 1 = OX(null(n_vars (p) ) ) ; 
for i=l to length(p) do 1 = 1+coef (p,i)X(expon(p,i)/k) od; 1 
Levi_branch2 (vec v; int i, j) = loc m = Levi_mat (i , j ) ; 
loc r = res_mat (m) ; loc p = r; 

for k = 1 to Lie_rank do p [k,Lie_rank] = od; p [i ,Lie_rank] = 1; 
for k = 1 to Lie_rank do p [k,Lie_rank-l] = od; p [j ,Lie_rank-l] = 1; 
loc q = (*p)*r; 

loc det_q = q[Lie_rank-l ,Lie_rank-l] *q[Lie_rank,Lie_rank] \ 
(2.12) -q[Lie_rank,Lie_rank-l] *q[Lie_rank-l ,Lie_rank] ; 

qq = q; 

for j = 1 to Lie_rank do qq[j ,Lie_rank] = -q[j ,Lie_rank] od; 

qq[Lie_rank-l ,Lie_rank] = 0; qq[Lie_rank,Lie_rank] = 1; 

for j = 1 to Lie_rank do qq[j ,Lie_rank-l] = -q[j ,Lie_rank-l] od; 

qq[Lie_rank-l ,Lie_rank-l] = 1; qq[Lie_rank,Lie_rank-l] = 0; 

loc s = null(Lie_rank,Lie_rank) ; 

for j = 1 to Lie_rank do s[j,j] = det_q od; 

s [Lie_rank-1 ,Lie_rank-l] = (q*qq) [Lie_rank,Lie_rank] ; 

s [Lie_rank,Lie_rank-l] = -(q*qq) [Lie_rank,Lie_rank-l] ; 

s [Lie_rank-l ,Lie_rank] = - (q*qq) [Lie_rank-1 ,Lie_rank] ; 

s [Lie_rank,Lie_rank] = (q*qq) [Lie_rank-1 ,Lie_rank-l] ; 

loc i_q = qq*s; 

loc b = branch(v, Cartan_type (m) , r) ; div(b*i_q* (*p) ,det_q) 

and to improve Levi_branch3 . lie is left as an exercise (which implicitly requires incorporating 
the formula for the inverse of a general 3x3 matrix). 



2D Case g simple, rank t = rankg and t is semisimple. 

This is the most delicate case: g is simple, 6 is an inner automorphism, and t is semisimple. Let us 
assume that 6 is ^-maximal. Then the Borel-de-Siebenthal structure theory [2] provides a simple 
root system ^ = {ipi, . . . , tp n } and a simple root j = ip r £ ^/ such that = (^ \ {7}) U {— /3 } is 
a simple root system for t, and 6 has order n r , where the maximal root /3 g = ^n^j- 

Let 5 denote the maximal rank subalgebra of g with simple root system ty s = (\& \ {7}). Let 
w g and w B denote, respectively, the longest elements of the Weyl groups W s and W s . We write 
S + (g,t) for the positive root system of q relative to t defined by 

Lemma 2.13 The transformation —w s preserves $ s and sends — (5 a into S + (g,t). 

Proof. In general, the longest element of the Weyl group sends the positive Weyl chamber to its 
negative, so — w s preserves <3? s . But — w s {— f3 s ) = w s ((3 g ) is obtained from f3 g = ^rriiipi by a series 
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of simple root reflections : £ i— > £ — j -0 with ^ ^ 7. Thus the coefficient of 7 in w s (/? g ) is the 
same as that in /3 S , which is n r > 0, so w s ((3 s ) £ £ + (g,t). □ 

We now indicate how the LiE program uses w s and (3 S to compute the restriction matrix res_wt, 
which it uses to calculate restrictions of representations of g to t. First, we use — w B to carry the 
simple root system = ^ s U {— (3 g } of t to another simple root system <3? := i 6 B U {w s ((3 s )}. The 
point is that <3? then consists of positive roots for g, all but one of them simple, by Lemma 12.131 
The LiE program assumes Bourbaki root order for both g and t. It permutes the roots of $ in a 
somewhat arbitrary way in order to do this when it computes the restriction matrix and applies it 
to branching of representations from g to t. We will try to do this in a way that involves minimal 
permutation. 

We start by computing w s within the LiE program. It is the long word for W s , but there LiE 
orders the roots incorrectly, so we use the slightly convoluted routine 

ws = reduce (long_word~r_reduce(long_word, [1,2, .. . ,r-l,r+l, . . . ,n-l,n] )) 

Here long_word is the longest element w s of the Weyl group W s , so the shortest element of the 
coset w s W e is r_reduce (long_word, [1,2, .. . ,r-l,r+l, . . . ,n-l,n]). Then we set up the new 
simple root system $ for t as the rows of a matrix RR in which w s (P s ) replaces 7 and the roots are 
re-ordered (minimally) to Bourbaki order. Then there are three ways to compute the restriction 
matrix res_wt. 

The first is to note that RR is the inverse of the matrix res_rt, so one can compute 
res_wt = i_Cartan(a) *rcs_rt*Cartan({) /det_Cartan(g) , 
which is (|2.5p . The second is just to use the LiE assignment res_wt = res_mat (RR) . And the third, 
which is in fact the way that LiE implements res_mat, is to set g as the default by setdef ault (g) 
(putting in the Lie type of g), initialize res_wt as a the n x n identity matrix, res_wt = id(n), 
and then fill it in by 

for = 1 to n do 

for i = 1 to n do 

res_wt[i,j] = Cartan(i_Cartan [i] ,RR[j] )/det_Cartan 

od 

od 

In the next few sections we will run through the various basic cases cases of (g, 6) and then 
reduce the general case to these basic cases. 

3 Cases: g is simple, t is (9— maximal and rankH < rankg 

Recall the automorphism 6 of g with t = g e . In this section we assume that t is ^-maximal, in 
other words that it is maximal among the proper ^-invariant subalgebras of g, and we apply the 
methods of Section I2BI 

Note that 9 is an outer automorphism of g because rankt < rankg. If some power 9 m ^ 1 is 
an inner automorphism then its fixed point set is ^-invariant and satisfies t ^ g e ™ ^ g. As 6 is 
^-maximal we conclude that every power 9 m 7^ 1 is an outer automorphism of g. All possibilities 
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are listed in [7J Theorem 5.10(3)]. There 9 has prime order p = 2 or p = 3. If p = 2 then G/if is 
one of the riemannian symmetric spaces 

SU{n)/SO(n), SU{2n)/Sp{n), SO(2p + 2 + 2q)/{SO{2p + 1) x SO{2q + 1)}, E 6 /F 4 , E & /Sp{A). 

If p = 3 then Cr/if is one of the nearly-kaehler spaces 

Spin(8)/G 2 and Spin(8)/SU(3). 

It is useful to note that either the Dynkin diagram of t is obtained by folding the diagram of g as 
in [IT] — all possibilities are listed in the tables of pp. 245, 247] — or 9 has form 9' o Ad(g) 
where g e is obtained by folding and g <E G e . For example Spin(8)/G2 is obtained by folding and 
Spin(8)/ SU(3) is derived from it as just described; and E^/F^ is obtained by folding and Eq/C^ 
is derived from it as just described. For details of the latter see |12l p. 291]. We now run through 
that list. 

3A Case G/K = SU "(2m)/ SO (2m). 

In order to find the restriction matrix used by the LiE program, we consider the Cartan subalgebras 
s = {diagjni, . . . , u m , —u m , . . . , — u{\ \ ui G V — 1 ^} of t and 
t = {diag{ui, . . . , u 2m } | Ui £ y/^lR, u\ H h u 2m = 0} of g. 

The simple roots of g are the ipi = Ei — £j+i for 1 ^ i < 2m, and the simple roots of 6 are the 
= — (for 1 ^ i < m) and <^ m = e m _i + e m . Thus the simple roots of g have restriction to 
s given by ipi\ s = <pi and ip m+ i\ s = (p m -i for 1 ^ i < m and ip m \ s = 2e m = (p m - <p m -i. Here is the 
relevant LiE routine branch_A_D . lie for branching from SU(2m) to SO(2m). It takes arguments 
(m,v), where m > is an integer and v is a vector of length 2m — 1 consisting of non-negative 
integers, and branches v. If m is already set in LiE then only the argument v is needed. 

# file branch_A_D . lie # 

# usage: branch_A_D(m, v) branches v from SU(2m) to SD(2m) # 

# and branch_A_D(v) does the same if m is already defined in LiE # 
branch_A_D(int m;vec v) = setdef ault(Lie_group(l,2*m - D); 
res_rt = null (2*m-l ,m) ; 

for i=l to m do res_rt[i,i] = 1 od; 

res_rt [m,m-l] = -1; res_rt[m,m] = 1; 

for i=l to m-1 do res_rt [m+i ,m-i] = 1 od; 

res_wt = i_Cartan*res_rt*Cartan(Lie_group(4,m) )/det_Cartan; 
answer = branch(v,Lie_group(4,m) ,res_wt) ; 
(3.1) printC'the branching of "+v+" from SU("+2*m+") to S0("+2*m+") is"); 

answer 

branch_A_D(vec v) = setdefault(Lie_group(l,2*m - D); 

res_rt = null (2*m-l ,m) ; 

for i=l to m do res_rt[i,i] = 1 od; 

res_rt [m,m-l] = -1; res_rt[m,m] = 1; 

for i=l to m-1 do res_rt [m+i ,m-i] = 1 od; 

res_wt = i_Cartan*res_rt*Cartan(Lie_group(4,m) )/det_Cartan; 
answer = branch(v,Lie_group(4,m) ,res_wt) ; 

printC'the branching of "+v+" from SU("+2*m+") to S0("+2*m+") is"); 
answer 
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Here is an example of its use: 



read branch_A_D . lie 

branch_A_D (8 ,[1,1,0,0,0,0,0,0,0,0,0,0,0,0,0]) 

the branching of [1,1,0,0,0,0,0,0,0,0,0,0,0,0,0] from SU(16) 
to S0(16) is 1X[1, 0,0, 0,0, 0,0,0] +1X [1 , 1 , , , , ,0 ,0] 

m=5 

branch_A_D( [1,2,0,0,0,0,0,0,0]) 

the branching of [1,2,0,0,0,0,0,0,0] from SU(10) to S0(10) is 
1X[1,0,0,0,0] +1X[1, 1,0,0,0] +1X[1,2,0,0,0] +1X[3,0,0,0,0] 

3B Case G/K = SU(2m + l)/SO(2m+ 1). 

In order to find the restriction matrix we consider the Cartan subalgebras 

5 = {diagjui, . . . , u m , 0, — u m , . . . , —ui} \ Ui G y/— IK} of t and 

t = {diag{ni, . . . , u 2m +i} | Ui G ^/^TM, ui H h u 2m +i = 0} of g. 

The simple roots of q are the ipi = £{ — £j+i for 1 ^ % ^ 2m, and the simple roots of 6 are the 
<ft = £ j — (f° r 1 = * < m ) an d (fm = s m (short simple root). Thus the simple roots of g have 
restriction to s given by ipi\ s = ipi and ip m+ i\ s = <~p m +i-i for 1 ^ i ^ m. Here is the relevant LiE 
routine branch_A_B . lie for branching from SU(2m+ 1) to SO(2m+ 1). It takes arguments (m, v) , 
where m > is an integer and v is a vector of length 2m consisting of non-negative integers, and 
branches v. If m is already set in LiE then only the argument v is needed. 

# file branch_A_B . lie # 

# usage: branch_A_B(m, v) branches v from SU(2m+l) to S0(2m+1) # 

# and branch_A_B (v) does the same if m is already defined in LiE # 
branch_A_B(int m;vec v) = setdef ault(Lie_group(l,2*m)) ; 

res_rt = null(2*m,m); 

for i=l to m do res_rt[i,i] = 1 od; 

for i=l to m do res_rt [m+i ,m+l-i] = 1 od; 

res_wt = i_Cartan*res_rt*Cartan(Lie_group(2,m))/det_Cartan; 
answer = branch(v,Lie_group(2,m) ,res_wt) ; 
printC'the branching of "+v+" from SU("+(2*m+l)+") to \ 
(3.2) S0("+(2*m+l)+") is"); 

answer 

branch_A_B(vec v) = setdef ault(Lie_group(l,2*m)) ; 

res_rt = null(2*m,m); 

for i=l to m do res_rt[i,i] = 1 od; 

for i=l to m do res_rt [m+i ,m+l-i] = 1 od; 

res_wt = i_Cartan*res_rt*Cartan(Lie_group(2,m))/det_Cartan; 

answer = branch(v,Lie_group(2,m) ,res_wt) ; 

printC'the branching of "+v+" from SU("+(2*m+l)+") to \ 

S0("+(2*m+l)+") is"); 
answer 
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3C Case G/K = SU(2m)/Sp(m). 



This case is quite similar to the case of SU {2m)/ SO {2m) above. We consider the Cartan subalge- 
bras 

s = {diagjui, . . . , u m , —u m , . . . , —ui} \ Ui £ \J — 1 M} of t and 

t = {diag{«i, . . . , U2m} | Ui £ u\ H V u 2m = 0} of q. 

The simple roots of q are the ipi = Si — £j+i for 1 ^ i < 2m, and the simple roots of 6 are the 
<Pi = £i~ (for 1 ^ i < m) and <p m = 2e m . Thus the simple roots of q have restriction to s given 
by ipi\ s = <Pi and ip m+ i\ s = <p m -i for 1 ^ i < m and ip m \s = 2e m = tp m ■ Here is the relevant LiE 
routine branch_A_C . lie for branching from SU{2m) to Sp(m). It takes arguments (m,v), where 
m > is an integer and v is a vector of length 2m — 1 consisting of non-negative integers, and 
branches v. If m is already set in LiE then only the argument v is needed. 

# file branch_A_C.lie # 

# usage: branch_A_C(m, v) branches v from SU(2m) to Sp(m) # 

# and branch_A_C(v) does the same if m is already defined in LiE # 
branch_A_C(int m;vec v) = setdef ault(Lie_group(l,2*m - 1)); 
res_rt = null (2*m-l ,m) ; 

for i=l to m-1 do res_rt[i,i] = 1 od; 
res_rt[m,m] = 1; 

for i=l to m-1 do res_rt [m+i ,m-i] = 1 od; 

res_wt = i_Cartan*res_rt*Cartan(Lie_group(3,m) )/det_Cartan; 
answer = branch(v,Lie_group(3,m) ,res_wt) ; 
(3.3) print ("the branching of "+v+" from SU("+2*m+") to Sp("+m+") is"); 

answer 

branch_A_C(vec v) = setdefault(Lie_group(l,2*m - 1)); 

res_rt = null (2*m-l ,m) ; 

for i=l to m-1 do res_rt[i,i] = 1 od; 

res_rt[m,m] = 1; 

for i=l to m-1 do res_rt [m+i ,m-i] = 1 od; 

res_wt = i_Cartan*res_rt*Cartan(Lie_group(3,m) )/det_Cartan; 
answer = branch(v,Lie_group(3,m) ,res_wt) ; 

print ("the branching of "+v+" from SU("+2*m+") to Sp("+m+") is"); 
answer 



3D Cases G/K = SO{2p + 2 + 2q)/{SO{2p + 1) x SO{2q+ 1)} (p,q not both 0). 

We use the Cartan subalgebras 

t = {u := diag{ui, . . . u p+q+ i, -u p+q+ i, m} \ Ui £ T 3 !!} and s = {u £ t | u p+q+1 = 0}. 

Now q has simple roots ipi = £i — £j+i for 1 ^ i ^ p + q and ip p+q+ i = e p+q + e p+q+ \. The 
subalgebra t has simple roots tpi = £\ — £j+i for 1 ^ i < p and for p + 1 ^ i < p + q, (p p = e p , and 
'■Pp+q = e p+q- Thus the simple roots of q have restriction to 5 given by ipi\ s = <pi for 1 ^ i < p and 
p <i <p + q, ipp\ s = (f p - Yli Pp+j, V'p+gls = <Pp+q, an d tpp +q+ i\ s = 'Pp+q- Here is the relevant LiE 
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routine branch_D_BB . lie for branching from SO(2p + 2q + 2) to SO(2p + 1) x SO(2q +1). It takes 
arguments (p , q, v) , where p, q ^ are integers not both zero and v is a vector of length p + q + 1 
consisting of non-negative integers, and branches v. If p and q are already set in LiE then only the 
argument v is needed. 

# file branch_D_BB.lie # 

# usage: branch_D_BB(p,q,v) branches v from S0(2p+2q+2) # 

# to S0(2p+l)xS0(2q+l) , and branch_D_BB(v) does the # 

# same if p and q are already defined in LiE # 
branch_D_BB (int p,q;vec v) = setdef ault (Lie_group(4,p+q+l) ) ; 
res_rt = null (p+q+1 ,p+q) ; 

for i=l to p-1 do res_rt[i,i] = 1 od; 

res_rt[p,p] = 1; for j=p+l to p+q do res_rt[p,j] = -1 od; 
for i=p+l to p+q-1 do res_rt[i,i] = 1 od; 
res_rt [p+q, p+q] = 1; res_rt [p+q+1 , p+q] = 1; 

res_wt = i_Cartan*res_rt *Cartan (Lie_group (2 , p) *Lie_group (2 , q) ) /4 ; 

answer = branch(v,Lie_group(2,p)*Lie_group(2,q) ,res_wt) ; 
. . print ("the branching of "+v+" from S0("+(2*p+2*q+2)+") to \ 

' ' S0("+(2*p+l)+")xS0("+(2*q+l)+") is"); 

answer 

branch_D_BB (vec v) = setdef ault (Lie_group(4, p+q+1) ) ; 

res_rt = null (p+q+1 , p+q) ; 

for i=l to p-1 do res_rt[i,i] = 1 od; 

res_rt[p,p] = 1; for j=p+l to p+q do res_rt[p,j] = -1 od; 
for i=p+l to p+q-1 do res_rt[i,i] = 1 od; 
res_rt [p+q, p+q] = 1; res_rt [p+q+1 , p+q] = 1; 

res_wt = i_Cartan*res_rt*Cartan (Lie_group (2 , p) *Lie_group (2 , q) ) /4 ; 
answer = branch(v,Lie_group(2,p)*Lie_group(2,q) ,res_wt) ; 
print ("the branching of "+v+" from S0("+(2*p+2*q+2)+") to \ 

S0("+(2*p+l)+")xS0("+(2*q+l)+") is") ; 
answer 



3E Case G/K = Spin(8)/G 2 . 

In this case g has simple root system {tpi, ip2, ip3, V^} numbered as in the Introduction, and t has 
simple root system {ipi, ^2} as follows. The Cartan subalgebra s of t is the subspace of the Cartan 
subalgebra t of g given by ipi = 1P3 = ipi- See [10]. The root restrictions are ipi\ s = = ^4 Is = ^1 
(the short simple root of 6) and ^2 Is = P2 (the long simple root of 6). Now the Lie routine for 
branching from Spin(8) to G2 is 

# file branch_D4_G2.1ie # 
branch_D4_G2(vec v) = setdef ault (D4) ; 
res_rt = [ [1 ,0] , [0 , 1] , [1 ,0] , [1 ,0] ] ; 
(3.5) res_wt = i_Cartan*res_rt*Cartan(G2) )/det_Cartan; 

answer = branch(v,G2,res_wt) ; 

printC'the branching of "+v+" from Spin(8) to G2 is"); 
answer 
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3F Cases G/K = Spin{8)/SU(3), G/K = E 6 /F± and G/K = E 6 /Sp(A). 

There the restriction matrices can be computed as in the case of Spin(8)/G2, or one can use the 
small database of maximal subalgebras built into LiE. That database is accessible by the commands 
resjnat (A2 ,D4) , resjnat (F4,E6) and resjnat (C4,E6) 

The corresponding LiE routines are 

# file branch_D4_A2.1ie # 
branch_D4_A2(vec v) = setdef ault (D4) ; 
answer = branch (v, A2, res_mat (A2, D4) ) ; 

printC'the branching of "+v+" from D4 to the triality A2 is"); 
answer 



# file branch_E6_F4 . lie # 
branch_E6_F4(vec v) = setdef ault (E6) ; 
answer = branch (v,F4,res_mat (F4,E6) ) ; 

printC'the branching of "+v+" from E6 to F4 is"); answer 

# file branch_E6_C4.1ie # 
branch_E6_C4(vec v) = setdef ault (E6) ; 
answer = branch (v, C4, res_mat (C4, E6) ) ; 

printC'the branching of "+v+" from E6 to Sp(4) is"); answer 



4 Cases: g is simple and t is the centralizer of a toral subalgebra 

These cases were covered in Section I2CI If t is the centralizer of a toral subalgebra of q it is the 
fixed point set of an automorphism. For if G is a connected Lie group with Lie algebra q and K is 
the analytic subgroup for t we can choose g E G such that the powers {g n \ n E Z} form a dense 
subgroup of the identity component of the center of K, and then t is the fixed point set of Ad(g). 

5 Cases: q is simple, 9 2 = 1 and rank£ = rankg 

In this section we apply the method of Section [2D] and run through the cases where n r = 2, i.e. the 
cases where G/K is a riemannian symmetric space. The other cases of equal rank will be considered 
in the next section. 

It turns out that, in the classical group symmetric space cases, we do not have to renumber the 
roots of fy s , while the renumbering is needed for most of the exceptional cases. 

5A Case (G, K) = (SO(2p + 2q+l), SO(2p) x SO(2q + 1)) where p ^ 2 and q ^ 0. 

Here w B reverses the order of tpi, . . . , ?/V-i but does not move tpi for i > r, so w s (f3 g ) attaches 
to the diagram of s at Vv-2- Now the simple roots of t in Bourbaki root order are given by 
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(3.6) 

and 
(3.7) 

and 
(3.8) 



{ipi, ■ ■ ■ ,ip r -i,w e (f3 g );ip r +i, . . . ,ip n }- Here is the LiE program branch_B_DB . lie that branches the 
representation of SO{2p + 2q + 1), specified by a vector v of length p + q, to SO{2p) x SO(2q + 1). 
Usage is branch J3_DB(p,q,v) , but if p and q are already set in LiE one can just use branch_B _DB (v) . 

#file branch_B_DB.lie # 
#usage: branch_B_DB(p,q,v) branches v from S0(2p+2q+l) to SD(2p)xSD(2q+l) # 
#and branch_B_DB (v) does the same if p and q are already defined in LiE # 

branch_B_DB(int p,q; vec v) = setdefault(Lie_group(2,p+q)) ; 

loc JJ = A1*A1; loc KK = Al; loc LL = Lie_group(4,p+q) ; 

u = null(p+q-l); for i = 1 to p+q-1 do u[i]=i od; 

if q > then for i = p to p+q-1 do u[i]=i+l od fi; 

ws = reduce (long_word~r_reduce (long_word,u) ) ; 

RR = id(p+q) ; RR[p] = W_rt_action(high_root ,ws) ; 

if p >= 3 then J J = Lie_group(4,p) fi; 

if q >= 2 then KK = Lie_group(2,q) fi; 

if q > then LL = JJ*KK fi; answer = branch (v, LL, res_mat (RR) ) ; 
print ("the branching of "+v+" from S0("+(2*p+2*q+l)+") " 
(5.1) "to S0("+2*p+")xS0("+(2*q+l)+") is "); 

answer 

branch_B_DB(vec v) = setdef ault (Lie_group(2 ,p+q) ) ; 

loc JJ = A1*A1; loc KK = Al; loc LL = Lie_group(4,p+q) ; 

u = null(p+q-l); for i = 1 to p+q-1 do u[i]=i od; 

if q > then for i = p to p+q-1 do u[i]=i+l od fi; 

ws = reduce (long_word~r_reduce (long_word,u) ) ; 

RR = id(p+q) ; RR[p] = W_rt_action(high_root ,ws) ; 

if p >= 3 then J J = Lie_group(4,p) fi; 

if q >= 2 then KK = Lie_group(2,q) fi; 

if q > then LL = JJ*KK fi; answer = branch (v, LL, res_mat (RR) ) ; 
print ("the branching of "+v+" from S0("+(2*p+2*q+l)+") " 

"to S0("+2*p+")xS0("+(2*q+l)+") is "); 
answer 



5B Case (G, K) = (SO(2p + 2q), SO(2p) x SO(2q)) where p, q ^ 2. 

Again w s reverses the order of ipi, . . . , Vv-i but does not move ipi for r < i ^ n — 2, so w s ((3 g ) 
attaches to the diagram of s at Vr-2 ( or doesn't attach, if r = 2). Now the simple roots of 
t in Bourbaki root order are {V>i, • • • , VV-i, w g (f3 g ); tp r +i, ■ ■ ■ > 4>n}- Here is the LiE program for 
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restriction of representations from SO(2p + 2q) to SO(2p) x SO{2q). 



#file branch_D_DD.lie # 
#usage: branch_D_DD(p,q,v) branches v from SD(2p+2q) to SD(2p)xS0(2q) # 
#and branch_D_DD(v) does the same if p and q are already defined in LiE # 
branch_D_DD(int p,q; vec v) = setdef ault (Lie_group(4,p+q) ) ; 
loc JJ = A1*A1; loc KK = A1*A1; 

u = null(p+q-l); for i = 1 to p+q-1 do u[i]=i od; 

for i = p to p+q-1 do u[i]=i+l od; 

ws = reduce (long_word~r_reduce (long_word,u) ) ; 

RR = id(p+q) ; RR[p] = W_rt_action(high_root ,ws) ; 

if p >= 3 then J J = Lie_group(4,p) fi; 

if q >= 3 then KK = Lie_group(4,q) fi; 

answer = branch (v, J J*KK,res_mat (RR) ) ; 

print ("the branching of "+v+" from S0("+(2*p+2*q)+") " 
(5.2) "to S0("+2*p+")xS0("+(2*q)+") is "); 

answer 

branch_D_DD(vec v) = setdef ault (Lie_group (4, p+q) ) ; 
loc JJ = A1*A1; loc KK = A1*A1; 

u = null(p+q-l); for i = 1 to p+q-1 do u[i]=i od; 

for i = p to p+q-1 do u[i]=i+l od; 

ws = reduce (long_word~r_reduce (long_word,u) ) ; 

RR = id (p+q) ; RR[p] = W_rt_action(high_root ,ws) ; 

if p >= 3 then J J = Lie_group(4,p) fi; 

if q >= 3 then KK = Lie_group(4,q) fi; 

answer = branch (v, J J*KK,res_mat (RR) ) ; 

print ("the branching of "+v+" from S0(" + (2*p+2*q) + ") 11 

"to S0("+2*p+")xS0("+(2*q)+") is "); 

answer 



5C Case (G, K) = (Sp(p + q), Sp(p) x Sp(q)) where p,q^.l. 

Again w s reverses the order of . . . ,ip r ~i but does not move ipi for i > r, so w s ({3 g ) attaches to 
the diagram of s at f/V-i (or doesn't attach, if r = 1). Now the simple roots of t in Bourbaki root 
order are {^i, . . . , Yv- 1, w s ((3 g ); ip r +i, ■ ■ ■ , V'n}- The LiE program for restriction of representations 
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from Sp(p + q) to Sp(p) x Sp(q) is 



#file branch_C_CC.lie # 
#usage: branch_C_CC(p,q,v) branches v from Sp(p+q) to Sp(p)xSp(q), # 
#branch_C_CC(v) does the same if p, q are already defined in LiE # 
branch_C_CC(int p,q; vec v) = setdefault(Lie_group(3,p+q)) ; 
loc JJ = Al; loc KK = Al; 

u = null(p+q-l); if p == 1 then for i = 1 to q do u[i]=i+l od fi; 
if p >= 2 then for i = 1 to p+q-1 do u[i]=i od fi; 
if p >= 2 then for i = p to p+q-1 do u[i]=i+l od fi; 
ws = reduce (long_word~r_reduce (long_word,u) ) ; 
RR = id(p+q) ; RR[p] = W_rt_action(high_root ,ws) ; 
if p >= 2 then JJ = Lie_group(3,p) fi; 
if q >= 2 then KK = Lie_group(3,q) fi; 
answer = branch (v, J J*KK,res_mat(RR) ) ; 
printC'the branching of "+v+" from Sp("+(p+q)+") " 
(5.3) "to Sp("+p+")xSp("+q+") is "); 

answer 

branch_C_CC(vec v) = setdef ault(Lie_group(3,p+q)) ; 
loc JJ = Al; loc KK = Al; 

u = null(p+q-l); if p == 1 then for i = 1 to q do u[i]=i+l od fi; 

if p >= 2 then for i = 1 to p+q-1 do u[i]=i od fi; 

if p >= 2 then for i = p to p+q-1 do u[i]=i+l od fi; 

ws = reduce (long_word~r_reduce (long_word,u) ) ; 

RR = id(p+q) ; RR[p] = W_rt_action(high_root ,ws) ; 

if p >= 2 then JJ = Lie_group(3,p) fi; 

if q >= 2 then KK = Lie_group(3,q) fi; 

answer = branch (v, J J*KK,res_mat(RR) ) ; 

printC'the branching of "+v+" from Sp("+(p+q)+") " 

"to Sp("+p+")xSp("+q+") is "); 

answer 



5D Case G 2 Mi^i- 

Now we run through the exceptional group cases. First suppose that G = G^- Then K = A±A\ 

with simple roots n=£n . Here /3 g = 3^i + 2ip 2 an d 7 = ip2, = {V"i}> an d w s(Pg) = Pg- Thus 
the LiE program (if one wants to bother with it in this case) is 

# file branch_G2_AlAl.lie # 

branch_G2_AlAl(vec v) = setdef ault (G2) ; 
, . ws = reduce (long_word~r_reduce (long_word, [1] )) ; 

^ ' ' RR = id(2); RR[2] = W_rt_action(high_root ,ws) ; 

answer = branch (v, A1A1, res_mat (RR) ) ; 

printC'the branching of "+v+" from G2 to A1A1 is"); answer 



20 



5E Cases F 4 / ' A X C Z and F 4 /5 4 . 



Ipl "02 ^3 "04 

Next suppose that G = F 4 . The simple root system is o — n-Y-n and the negative of the 

maximal root /3 S = 2ipi + 3^2 + 4?/>3 + 2"0 4 attaches at ■i/'i- Thus there are two possibilities: y = ifti 
and i'T = A1C3, or 7 = -0 4 and K = B4. The corresponding LiE programs are given by 

# file branch_F4_AlC3.1ie # 

branch_F4_AlC3(vec v) = setdef ault (F4) ; 

ws = reduce (long_word~r_reduce(long_word, [2,3,4] )) ; 

(5.5) RR = id(4); RR[1] = W_rt_action(high_root ,ws) ; 
RR[2] = id (4) [4] ; RR[4] = id (4) [2] ; 

answer = branch (v, A1C3, res_mat (RR) ) ; 

printO'the branching of "+v+" from F4 to A1C3 is"); answer 

and 

# file branch_F4_B4.1ie # 

branch_F4_B4(vec v) = setdef ault (F4) ; 

ws = reduce (long_word~r_reduce (long_word, [1 , 2 ,3] ) ) ; 

(5.6) RR = id(4); RR[1] = W_rt_action(high_root ,ws) ; 

RR[2] = id(4)[l]; RR[3] = id(4) [2] ; RR[4] = id(4) [3] ; 
answer = branch (v, B4, res_mat (RR) ) ; 

printO'the branching of "+v+" from F4 to B4 is"); answer 



5F Cases E % jA x A^ and E % jA h A x . 

Now suppose that G = Eq. Then there are three simple roots of coefficient 2 in the maximal 
root. All of them differ by automorphisms of the extended Dynkin diagram, so the corresponding 
subalgebras t differ by an automorphism, but in [5] we will need to distinguish between them , so 

i>\ "03 V>4 V>5 1p<5 

we treat them separately. The simple root system is J ^° and the negative of the 

maximal root (3 g = tpi + 2-02 + 2^3 + 3^4 + 2^5 + ipQ attaches at tp2- There are three equivalent 
possibilities: (i) 7 = ips and K = A\A^, (ii) 7 = V's an d K = A5A1, and (hi) 7 = ip2 and K = A$A\. 
For the first, the LiE routine is 

# file branch_E6_AlA5.1ie # 
branch_E6_AlA5(vec v) = setdef ault (E6) ; 

ws = reduce (long_word~r_reduce (long_word, [1,2,4,5,6] )) ; 

(5.7) RR = id(6); RR[6] = W_rt_action(high_root ,ws) ; 

RR[3] = id(6)[4]; RR[4] = id(6) [5] ; RR[5] = id(6) [6] ; 
answer = branch(v,AlA5,res_mat(RR)) ; 

printO'the branching of "+v+" from E6 to A1A5 is"); answer 
The second is quite similar, 

# file branch_E6_A5Al.lie # 
branch_E6_A5Al (vec v) = setdef ault (E6) ; 

ws = reduce (long_word~r_reduce (long_word, [1,2,3,4,6] )) ; 

(5.8) RR = id(6); RR[1] = W_rt_action(high_root ,ws) ; 
RR[2] = id(6)[l]; RR[5] = id(6) [2] ; 

answer = branch(v, A5A1 ,res_mat (RR) ) ; 

printO'the branching of "+v+" from E6 to A5A1 is"); answer 
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and the third is a bit different, 

# file branch_E6_A5Ala.lie # 
branch_E6_A5Ala(vec v) = setdef ault (E6) ; 
ws = reduce (long_word~r_reduce (long_word, [1,3,4,5,6] )) ; 
(5.9) RR = id(6); RR[6] = W_rt_action(high_root ,ws) ; 

RR[2] = id(6)[3]; RR[3] = id(6) [4] ; RR[4] = id(6) [5] ; RR[5] = id(6) [6] ; 
answer = branch(v, A5A1 ,res_mat (RR) ) ; 

printO'the branching of "+v+" from E6 to A5A1 is"); answer 



5G Cases E 7 /A x Dq, E 7 /D 6 A 1 and E 7 /A 7 . 

Next, let G = E 7 . Then there are three simple roots of coefficient 2 in the maximal root. Two of 
them differ by an automorphism of the extended Dynkin diagram, so the corresponding subalgebras 
6 differ by an automorphism, but we will need to distinguish between them in [5]. So, as in some 

■01 V3 V>4 4>5 ^6 1p7 



(5.10) 



of the Eq cases above, we treat them separately. The simple root system is I ip 2 

and the negative of the maximal root (5 S = 2ipi + 2ip2 + 3^3 + 4^4 + 3ip$ + 2-^6 + ^7 attaches at ipi . 
The possibilities are (i) j = ipi and K = A\D§, (ii) 7 = ipQ and K = DqAi, and (iii) 7 = 1P2 an d 
K = A 7 . For first of these the LiE routine is 

# file branch_E7_AlD6.1ie # 
branch_E7_AlD6 (vec v) = setdef ault (E7) ; 

ws = reduce (long_word~r_reduce (long_word, [2,3,4,5,6,7] )) ; 
RR = id(7); RR[1] = W_rt_action(high_root ,ws) ; 
RR[2] = id(7)[7]; RR[3] = id(7) [6] ; RR[4] = id(7) [5] ; 
RR[5] = id(7)[4]; RR[6] = id(7) [3] ; RR[7] = id(7) [2] ; 
answer = branch(v,AlD6,res_mat(RR)) ; 

printO'the branching of "+v+" from E7 to A1D6 is"); answer 

For obvious reasons the second is similar 

# file branch_E7_D6Al.lie # 
branch_E7_D6Al (vec v) = setdef ault (E7) ; 

ws = reduce (long_word~r_reduce (long_word, [1,2,3,4,5,7] )) ; 
(5.11) RR = id(7); RR[1] = W_rt_action(high_root ,ws) ; 

RR[2] = id(7)[l]; RR[6] = id(7) [2] ; 
answer = branch (v, D6A1 ,res_mat (RR) ) ; 

printO'the branching of "+v+" from E7 to D6A1 is"); answer 

and the third is a bit different 

# file branch_E7_A7 . lie # 
branch_E7_A7(vec v) = setdef ault (E7) ; 

ws = reduce (long_word~r_reduce (long_word, [1,3,4,5,6,7] )) ; 
, , RR = id(7) ; RR[7] = W_rt_action(high_root ,ws) ; 

^ > RR[2] = id(7)[3]; RR[3] = id(7) [4] ; RR[4] = id(7) [5] ; 

RR[5] = id(7)[6]; RR[6] = id(7) [7] ; 

answer = branch (v, A7, res_mat (RR) ) ; 

printO'the branching of "+v+" from E7 to A7 is"); answer 
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5H Cases E 8 /D 8 and E 8 /E 7 A 1 

tpl 1p3 1p4 1p5 V>6 1p7 1p8 

Finally, suppose that G = E 8 . The simple root system is \' l p2 an< ^ the negative 

of the maximal root f3 g = 2ipi + 3ip 2 + 4^3 + 6^4 + 5^5 + 4^6 + 3-07 + 2-i/>g attaches at ip 8 ■ The 
possibilities are (i) 7 = tpi and K = D 8 and (ii) 7 = ip 8 an d K = E 7 A\. In the first case the LiE 
routine is 

# file branch_E8_D8.1ie # 
branch_E8_D8(vec v) = setdef ault (E8) ; 

ws = reduce (long_word~r_reduce(long_word, [2,3,4,5,6,7,8] )) ; 
RR = id(8); RR[1] = W_rt_action(high_root ,ws) ; 
RR[2] = id(8) [8] ; RR[3] = id(8) [7] ; RR[4] = id(8) [6] ; 
RR[6] = id(8) [4] ; RR[7] = id(8) [3] ; RR[8] = id(8) [2] ; 
answer = branch ( v, D8, res_mat (RR) ) ; 

printO'the branching of "+v+" from E8 to D8 is"); answer 

and in the second it is 

# file branch_E8_E7Al.lie # 
branch_E8_E7Al(vec v) = setdef ault (E8) ; 

ws = reduce (long_word~r_reduce (long_word, [1,2,3,4,5,6,7] )) ; 
\ ' ' RR = id(8); RR[8] = W_rt_action(high_root ,ws) ; 

answer = branch (v, E7A1 ,res_mat (RR) ) ; 

printO'the branching of "+v+" from E8 to E7A1 is"); answer 



(5.13) 



Taking into account the results of Sections [3] and [H now for every compact connected, simply 
connected symmetric space G/K, G simple, we have shown how to compute the restriction to K of 
any irreducible finite dimensional representation of G. The irreducible compact connected, simply 
connected symmetric spaces G/K with q not simple are the are the simply connected simple Lie 
group manifolds, which were covered in Section [2AI 



6 Cases: q is simple, 9 3 = 1 or # 5 = 1, and rank£ = rankg 

The maximal connected subgroups of maximal rank in a compact connected Lie group were de- 
scribed by A. Borel and J. de Siebenthal [2j. Most of them are symmetric subgroups, and their 
classification can be used in the classification of symmetric spaces |12| . The ones that are sym- 
metric were considered in Section [5j The others correspond to the simple roots 7 whose coeffi- 
cient in the maximal root /3 g is an odd prime, necessarily 3 or 5. The ones for prime 3 are given 
by (G, K) = (G 2 , A 2 ), (F 4 , A 2 A 2 ), (E 6 ,A 2 A 2 A 2 ), (E 7 , A 2 A 5 ), (E 7 , A 5 A 2 ), (E 8 , A 8 ), and (E 8 , E 6 A 2 ). 
The one for prime 5 is given by (G,K) = (E 8 , A4A4). In all cases a simple root system for t is 
given by ^ = (\& g \ {7}) U {— /3 }, so we can use the methods of Section |2C1 as in Section [5j 

In all of these cases one can rely on a LiE database to produce a restriction matrix res_mat (Y , X) . 
However the applications in [S] require that we keep track of which root of Y comes from which root 
of X, and LiE scrambles the root order, so we generally have to do this by hand. In each case we 
indicate which elements of the simple root system <I> = {ipi, ■ ■ ■ , <p n } of 6 come from which elements 
of \ {7}) U {w B ((3 s )}, where we try to use the least complicated correspondence. 
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6A Case G/K = G 2 /A 2 . 



Here 7 = tpi, t is of type A 2 = su(3), and $ = {(fi, ^2} where if \ comes from w s (P s ) and if2 comes 
from ip2- This is indicated in terms of the Dynkin diagrams, by 



The corresponding LiE routine is 



a e i > o ^o. o o 

7 -h «>«(&) 



# file branch_G2_A2.1ie # 

branch_G2_A2(vec v) = setdef ault (G2) ; 
, . ws = reduce(long_word~r_reduce(long_word, [2] )) ; 

^ ' ' RR = id(2); RR[1] = W_rt_action(high_root ,ws) ; 

answer = branch (v, A2, res_mat (RR) ) ; 

printO'the branching of "+v+" from G2 to A2 is "); answer 



6B Case G/K = F A /A 2 A 2 . 

Here 7 = V>2, % is of type A2A2, and = {<fi,<P2,f3,f4} where ip\ comes from ipi, ip2 comes from 
w s ({3 g ), (f3 comes from ^3, and <p4 comes from ^4. This is indicated in terms of Dynkin diagrams 
by 

l/>l V>2 4>3 V>4 tfl tp2 ipz <fi4 
O O Q ) D O O O ) O O 

-Pt 7 10, (ft) 

The corresponding LiE routine is 

# file branch_F4_A2A2.1ie # 

branch_F4_A2A2(vec v) = setdef ault (F4) ; 
, . ws = reduce (long_word~r_reduce (long_word, [1 , 3,4] )) ; 

^ ' ' RR = id(4); RR[2] = W_rt_action(high_root ,ws) ; 

answer = branch (v, A2A2, res_mat (RR) ) ; 

printO'the branching of "+v+" from F4 to A2A2 is "); answer 



6C Case G/K = E 6 /A 2 A 2 A2. 



Here 7 = ip4, t is of type A2A2A2, and <E> = {(fi, <f2, <f3, <^4, <Pg} where <pi comes from ipi, if2 
comes from ^3, (p% comes from ip2, comes from w s ((3 s ), (p§ comes from ^5, and (p$ comes from 
ipQ. This is indicated in terms of Dynkin diagrams by 



'01 

o- 



V>3 



V>4 
— O- 



1p5 



V>6 
— O 



7 



^20 



W s (/9 fl ) 
— O V4 O 



¥>6 



¥>6 
^3 



¥>3 6 



6-/3 B 
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The corresponding LiE routine is 

# file branch_E6_A2A2A2.1ie # 
branch_E6_A2A2A2(vec v) = setdef ault (E6) ; 
ws = reduce (long_word~r_reduce (long_word, [1,2,3,5,6] )) ; 
(6.3) RR = id(6); RR[4] = W_rt_action(high_root ,ws) ; 

RR[2] = id(6)[3]; RR[3] = id(6) [2] ; 
answer = branch (v, A2A2A2, res_mat (RR) ) ; 

printO'the branching of "+v+" from E6 to A2A2A2 is"); answer 



6D Case G/K = E 7 /A 2 A 5 . 

Here 7 = ips, £ is of type A2A5, and <E> = {<pi, (f2, <^3, <^4, <^5, <f6, ^7} where comes from ipi, ip2 
comes from w s ({3 g ), ips comes from ip2, <p4 comes from ^4, ip$ comes from ^5, V?6 comes from tpQ, 
and ip-j comes from ip-j. This is indicated in terms of Dynkin diagrams by 

Ipl i>3 i>4 i>5 i>6 4>7 Vl V2 <f4 ¥>5 ¥>6 <P7 
O O O O O O O O O O O O O 

-h 7 ^ ] 

l/<2 O ¥>3 O 

The corresponding LiE routine is 

# file branch_E7_A2A5.1ie # 
branch_E7_A2A5(vec v) = setdef ault (E7) ; 

ws = reduce (long_word~r_reduce (long_word, [1,2,4,5,6,7] )) ; 
(6.4) RR = id(7); RR[2] = W_rt_action(high_root ,ws) ; RR[3] = id(7) [2] ; 

res = res_mat (RR) ; 

answer = branch ( v, A2A5, res_mat (RR) ) ; 

printC'the branching of "+v+" from E7 to A2A5 is"); answer 



6E Case G/K = E 7 /A 5 A 2 . 



Here 7 = ^5, t is of type A5A2, and = {<pi, (f2, ^3; ^4> ¥>6j where </?i comes from ^1, 992 
comes from ^3, 923 comes from ^4, 924 comes from 1^2, 925 comes from w s (/3 g ), (p§ comes from ^/>6> 
and (f7 comes from ip-j. This is indicated in terms of Dynkin diagrams by 



■•ci- 



v's 



■ij: 4 



V>5 

7 



V>6 



V>7 



o- 



<P2 
—O- 



V>3 



■02 O 



<P6 
O- 



<P7 
— O 



■-P4 
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The corresponding LiE routine is 

# file branch_E7_A5A2 . lie # 
branch_E7_A5A2 (vec v) = setdef ault (E7) ; 

ws = reduce (long_word~r_reduce (long_word, [1,2,3,4,6,7] )) ; 
(6.5) RR = id(7); RR[5] = W_rt_action(high_root ,ws) ; 

RR[2] = id(7)[3]; RR[3] = id(7) [4] ; RR[4] = id(7) [2] ; 
answer = branch(v,A5A2,res_mat(RR)) ; 

printC'the branching of "+v+" from E7 to A5A2 is"); answer 



6F Case G/K = E 8 /A 8 . 

Here 7 = ip2, 6 is of type A$, and & = {ifi, f2, f3, if 4, <P5, Pe, ^7, Vs} where f\ comes from w s ((3 s ), 
if 2 comes from ipi, if 3 comes from ^3, fi comes from ^4, ip$ comes from ^5, fQ comes from tpQ, ipj 
comes from ^7, and comes from ip 8 . This is indicated in terms of Dynkin diagrams by 

V>1 1p3 1p4 1p5 i>6 1p7 4>8 <P2 f.i V>4 lp5 <P6 ipi fS 
O O O O O O O O O O O O O O O O 

-/?„ 10,09,) 

ip 2 o 7 

The corresponding LiE routine is 

# file branch_E8_A8.1ie # 

branch_E8_A8(vec v) = setdef ault (E8) ; 
. . ws = reduce (long_word~r_reduce (long_word, [1 ,3 ,4, 5 ,6,7,8] ) ) ; 

' ' ' RR = id(8); RR[1] = W_rt_action(high_root ,ws) ; RR[2] = id(8) [1] ; 

answer = branch (v, A8, res_mat (RR) ) ; 

printC'the branching of "+v+" from E8 to A8 is"); answer 



6G Case G/K = E 8 /E 6 A 2 . 

Here 7 = ipj, t is of type EqA2, and $ = {ipi, if 2, V>3, y>4, <P5, fe, <P7, Vs} where ipi comes from ipi, 
f2 comes from ip2, ¥3 comes from ips, if 4 comes from ^4, ip^ comes from ^5, ipg comes from ip§, if-j 
comes from w s ((3g), and <f 8 comes from ip 8 . This is indicated in terms of Dynkin diagrams by 

Ipi 1p:i 4>A V>5 1pd 1p7 i>8 Vl <P3 ¥>4 <P5 <P6 f7 <P8 

o o o o o o o o o o o o o o o 

7 -P a ~* \ w s (M 

1p2 O ¥2 O 

The corresponding LiE routine is 

# file branch_E8_E6A2.1ie # 

branch_E8_E6A2(vec v) = setdef ault (E8) ; 
, ws = reduce (long_word~r_reduce (long_word, [1,2,3,4,5,6,8] )) ; 

^ ' ' RR = id(8); RR[7] = W_rt_action(high_root ,ws) ; 

answer = branch (v, E6A2, res_mat (RR) ) ; 

printC'the branching of "+v+" from E8 to E6A2 is"); answer 
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6H Case G/K = E 8 /A±A 4 . 



Here 7 = fa, t is of type A4A4, and <3? = {</?i, <^2> V?3j fit ¥>5> V6) y7> Vs} where (pi comes from fa, 
if 2 comes from fa, cp^ comes from fa, (^4 comes from fa, f§ comes from w s {(3 s ), (p§ comes from 
ipQ, (pi comes from fa, and </?s comes from This is indicated in terms of Dynkin diagrams by 

1pl 1p3 l/>4 i>5 1p& Ipl i>& 'Pi V2 ¥3 fh <PB <P7 <P8 
O O O O O O O O O O O O O O O 

1p2 O <P4 

The corresponding LiE routine is 

# file branch_E8_A4A4.1ie # 
branch_E8_A4A4(vec v) = setdef ault (E8) ; 

ws = reduce (long_word~r_reduce (long_word, [1,2,3,4,6,7,8] )) ; 
(6.8) RR = id(8); RR[5] = W_rt_action(high_root ,ws) ; 

RR[2] = id(8)[3]; RR[3] = id(8) [4] ; RR[4] = id(8) [2] ; 
answer = branch (v, A4A4,res_mat (RR) ) ; 

printO'the branching of "+v+" from E8 to A4A4 is"); answer 
This completes our branching project as described in Section 2. 
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